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ABSTRACT 

Maintaining the didactic objective of the CEINCI-LAB Computer system, 
this article presents new functions that facilitate the spatial seismic analysis of an 
irregular structure in plan and elevation, generates the three-dimensional drawing 
of the structure, the data entry of the sections of its elements in structures formed 
only by beams and columns with unique sections. 

Performing the calculation of the stiffness matrix in floor coordinates, is 
one of the most important and complicated aspects of seismic analysis and that is 
why the flat and spatial analysis of an irregular structure will be exemplified. 

keywords: New functions of CEINCI-LAB. Spectral method. 3D Irregular 
Structures. 

NUEVAS FUNCIONES DEL SISTEMA CEINCI-LAB 
PARA ANÁLISIS SÍSMICO ESPACIAL DE 
ESTRUCTURAS IRREGULARES 

RESUMEN 

Manteniendo el objetivo didáctico dei sistema de computación CEINCI- 
LAB, en este artículo se presentan nuevas funciones que facilitan el análisis 
sísmico espacial de una estructura irregular en planta y elevación, genera el dibujo 
en tres dimensiones de la estructura, la entrada de datos de las secciones de sus 
elementos en estructuras conformadas solo por vigas y columnas con secciones 
únicas. 


Realizar el cálculo de la matriz de rigidez en coordenadas de piso, es uno 
de los aspectos más importantes y complicados dei análisis sísmico y es por ello 
por lo que se ejemplificará el análisis plano y espacial de una estructura irregular. 

Palabras Claves: Nuevas funciones de CEINCI-LAB. Método espectral. 
Estructuras espaciales irregulares. 
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1. INTRODUCCIÓN 

El usuário dei sistema de computación CEINCI-LAB debe conocer cómo 
se resuelve una estructura ante cargas estáticas o la forma de realizar ei análisis 
sísmico de un edificio en tres dimensiones. 

Para cumplir con ei objetivo dei análisis sísmico, en este artículo, se 
presentan nuevas funciones que permiten graficar y analizar sísmicamente una 
estructura irregular en ei plano y en tres dimensiones, la entrada de datos es 
sencilla y por ello los usuários se van a motivar en realizar ei análisis sísmico 
empleando las rutinas de CEINCI-LAB. 

CEINCI-LAB es un programa didáctico cuyo espíritu es ayudar a los 
usuários a entender la secuencia de cálculo dei análisis estructural y con ello 
apoyar a afianzar sus conocimientos con la realización de ejercicios; sin embargo, 
con las nuevas contribuciones se facilita ei ingreso de la geometria de la 
estructura. 

En artículos anteriores los autores han presentado funciones para ei 
análisis sísmico de estructuras en ei plano y espaciales regulares. Este artículo se 
convierte en un complemento a estos trabajos, pero para estructuras espaciales 
irregulares con pórticos conformados por columnas y vigas de secciones únicas. 

2. FUNCIONES DE CEINCI-LAB 

Es recomendable que ei lector conozca las funciones básicas que se 
encuentran en Aguiar (2014) y las nuevas funciones que se presentaron en Aguiar 
et ai (2019), estas últimas muy rápidamente son: 

• [nv,np,nudt,nudcol,nudvg,nod,nr]=geometria_nudo_viga(sv,sp) 


El programa geometria_nudo_viga facilita notablemente ei ingreso de 
datos de la geometria de un pórtico regular; los datos de ingreso son: sv que es 
un vector que contiene las luces de cada uno de los vanos; sp es otro vector con 
la altura de cada uno de los pisos. El programa reporte: nv número de vanos; np 
número de pisos; nudt es ei número de elementos totales, se recuerda que en 
cada viga se considera un nudo en la mitad; nudcol es ei número de columnas; 
nudvg es ei número de vigas; nod es ei número de nudos; nr es ei número de 
nudos restringidos. 

• [X,Y]=glinea_portico2(nv,np,sv,sp,nod,nr) 


Este programa es similar a la función glinea_portico (ver funciones 
básicas de CEINCI-LAB, en Aguiar, 2014), pero la entrada de datos es diferente. 
Determinar los dos vectores X, Y con parte de los resultados que reporta ei 
programa anterior. Solo sirve para pórticos regulares, considerando nudo en la 
mitad de las vigas. 
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• [NI,NJ]=gn_portico2(nr, nv, nudt, nudcol, nudvg) 


Programa similar a gn j oortico que determina los vectores con las 
coordenadas dei nudo inicial y final denominados NI,NJ pero en esta ocasión los 
datos de entrada son diferentes. 

• [X, Y,Z,nvgY,nportz,nodA ]=glinea_portico3D(CoordY,nr,np,nv,sp,nod,X, Y) 

Este programa genera las coordenadas de los nudos: X, Y, Z de la 
estructura espacial. El usuário primero indica la geometria de los Pórticos en 
sentido X; en efecto, para cada pórtico da las luces entre vanos (sv) y la altura de 
piso (sp). Posteriormente se requiere especificar los pórticos en las coordenadas Y 
para ello se emplea el vector denominado CoordY, se recomienda imponerse 
unos ejes de referencia X, Y en la parte inferior izquierda de la estructura. 

• [NI, NJ]=gn_portico3D(nudcol, nudvg, nvgY,nportz,nod, nudt, nodA, NI, NJ) 


Programa para generar el Nudo inicial y final de los elementos 
denominados NI, NJ pero de la estructura en tres dimensiones. Los datos de 
ingreso son: nudcol número de columnas; nudvg número de vigas; nvgY Número 
de vigas en Y; nportz Número de pórticos en Y; nod número de nudos; nudt 
número de elementos; nodA Número de nudos adicionales (nudo medio de las 
vigas en Y); NI, NJ Vectores con los nudos iniciales y finales generados, para un 
pórtico en el plano XZ. 

• dibujo3D(X, Y,Z,NI,NJ) 


Programa para dibujar una estructura tridimensional, este programa ayuda 
a visualizar en 3 dimensiones la estructura. Los datos para ingresar son: X Vector 
que contiene coordenadas en X; Y Vector con las coordenadas en Y, Z Vector que 
contiene coordenadas en Z; NI Vector con los nudos iniciales de los elementos; NJ 
Vector con los nudos finales de los elementos. 


3. NUEVAS FUNCIONES DE CEINCI-LAB 

En esta sección dei artículo se presentan las nuevas funciones para 
graficar la estructura irregular en dos y tres dimensiones, para el ingreso de las 
secciones de los elementos y para el análisis sísmico espacial considerando un 
modelo de piso rígido con tres grados de libertad por planta. 

• [VecComb]=Combinar_ VectoresfNum Vect, VI, V2, V3, V4, V5, V6, V7, V8, V9, V10) 


Este programa combina vectores en una matriz VecComb, el propósito es 
mantener la información referente a la geometria de la estructura de forma 
compacta. Los datos son: NumVect el número de vectores a combinar; VI, V2,..., 
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V10 son los vectores con la información a combinar (no requieren ser de la misma 
dimensión); en caso de requerir combinar menos de 10 vectores se debe colocar 
cero en lugar de Vi. 

• [nr,nod,np,nELEM,X, Y,NI,NJ]=geom_2D_lrreg(Num_Comb,svF,spF,Xo 
F,YoF); 

Este programa permite determinar la geometria de un pórtico plano (2D), 
calcula parâmetros necesarios para el análisis sísmico como nr número de 
restricciones; nod número de nudos; np número de pisos, nELEM cantidad de 
elementos; también se presenta las coordenadas de los nudos: X, Y de la 
estructura en el plano en los vectores X, Y respectivamente; NI Vector con los 
nudos iniciales de los elementos; NJ Vector con los nudos finales de los 
elementos. Se requiere de datos la siguiente información: Num_Comb tiene que 
ver con el número de partes que componen a la estructura; svF matriz que 
contiene las luces de cada uno de los vanos de los pórticos correspondientes a las 
“partes que componen la estructura” y de forma similar la matriz spF pero 
referente a las alturas de piso; XoF, YoF que permiten desfasar en sentido X o Y a 
las partes que componen a la estructura respecto a los ejes de referencia X, Y dei 
origen de coordenadas impuestos. 

• [nPortX,nPortY,X,Y,Z,NI,NJ]=geom_3D_lrreg(Num_Comb,svF,spF,Coo 
rdYF,XoF,ZoF) 


Programa para identificar el número de pórticos que componen la 
estructura espacial en los sentidos “XZ’ nPortX y “YZ’ nPortY, también se 
presenta las coordenadas de los nudos: X, Y, Z de la estructura espacial en los 
vectores X, Y, Z respectivamente; NI Vector con los nudos iniciales de los 
elementos; NJ Vector con los nudos finales de los elementos. Los datos 
requeridos son similares a los detallados en la rutina anterior geom_2D_lrreg pero 
con la información de la estructura espacial, más adelante en el ejemplo se podrá 
comprender de mejor manera a lo que se refiere con “partes que componen a la 
estructura”. La matriz CoordYF con las coordenadas de los pórticos en sentido Y 
de cada parte de la estructura es lo único que difiere en el caso espacial. 

• dibujo_2Di(X, Y,Z,NI,NJ,nelem,dibuj_plan) 


Programa para dibujar pórticos en 2D de una estructura tridimensional y 
genera la vista en planta de esta. Los datos para ingresar son: X Vector que 
contiene coordenadas en X; Y Vector con las coordenadas en Y, Z Vector que 
contiene coordenadas en Z; NI Vector con los nudos iniciales de los elementos; NJ 
Vector con los nudos finales de los elementos; nelem número de elementos de la 
estructura principal (columnas y vigas), en estructuras con elementos diagonales 
la rutina representa a los elementos diagonales de color rojo y a la estructura 
principal de color azul; dibuj_plan es una matriz de 2 filas, se identifica a la 
primera columna con cero o uno de acuerdo ai sentido dei pórtico que se desea 
representar, donde cero corresponde ai sentido XZ y uno para el sentido YZ, a 
continuación de estos números se escribe una lista con el número de pórtico que 
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se desea dibujar en cada sentido. Es importante indicar que las dos filas de la 
matriz dibuj_plan debe contener la misma cantidad de columnas, en el caso de 
querer representar más pórticos en un sentido que en otro se debe completar con 
ceros los espacios para completar las mismas columnas en ambas filas. 

• dibujoNudElem(X, Y,NI,NJ,CG) 


Esta función permite graficar pórticos planos en cada iteración de análisis 
para los sentidos XZ o YZ y presenta de color rojo los nudos con la identificación 
de cada nudo (nombre o etiqueta asignado en cada análisis, siempre empieza a 
enumerar en uno) además de color azul identifica a cada elemento dei pórtico. Los 
datos requeridos son: X Vector que contiene coordenadas en X; Y Vector con las 
coordenadas en Y (X representa las coordenadas e los nudos en sentido de 
análisis, Y representan las coordenadas en elevación de los nudos); NI Vector con 
los nudos iniciales de los elementos; NJ Vector con los nudos finales de los 
elementos; CG grados de libertad asignados a los nudos, toda la información debe 
corresponder al pórtico de análisis. 

• [nglXi,ELEMXi,Lxi,senoxi,cosenoxi,VCXi,Xxi,Zxi,NIXi,NJXi,CGXi]=Dat 
os_KLx_3D_lrreg(i,nPortY,ELEMX,X,Y,Z,NI,NJ); 


Este programa permite generar los datos para calcular la rigidez de 
pórticos planos en sentido XZ en la iteración para pórticos 3D. Los datos a 
ingresar son: i Indica el pórtico de análisis (Ej. i=1, Pórtico 1 y si i=n, Pórtico n); 
nPortY Número de pórticos en X; ELEMX Matriz que contiene la base y la altura 
de cada elemento o área e inércia si es acero; X coordenadas "X" de los nudos; Y 
coordenadas "Y" de los nudos, Z coordenadas en “Z’ de los nudos; NI, NJ 
Vectores con los nudos iniciales y finales generados. 

Lo que nos reporta el programa es lo siguiente: ngIXi Número de grados 
de libertad en el pórtico de análisis; ELEMXi, Matriz que contiene la base y la 
altura de cada elemento o área e inércia si es acero, pero dei pórtico de análisis; 
Lxi Vector que contiene la longitud de los elementos; senoxi Vector que contiene 
los senos de los elementos; cosenoxi Vector que contiene los cosenos de los 
elementos; VCXi Matriz que contiene los vectores de colocación de elementos, Xxi 
coordenadas "X" de los nudos; Zxi coordenadas "Z' de los nudos, Nlxi, NJxi 
Vectores con los nudos iniciales y finales generados y CGXi los grados de libertad 
asignados a los nudos, nuevamente todo esto dei pórtico de análisis en sentido X. 

Esta rutina identifica los nudos y elementos que corresponden al pórtico de 
análisis en sentido XZ y reordena la información necesaria específicamente para el 
análisis de dicho pórtico en nuevos vectores o matrices identificados como “..xi”. 

La información referente a la sección de los elementos en la matriz 
ELEMX solo corresponde a columnas y vigas, es decir una matriz cuadrada de 
2x2 donde la primera fila corresponde a la base y altura de columnas y la segunda 
fila lo propio para vigas; en caso de ser de acero se debe ingresar área e inércia 
de los elementos de las columnas y de vigas. En la rutina se identifica los 
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elementos verticales como columnas a los cuales se le asigna la propiedad de 
columnas, luego se analiza que los elementos horizontales corresponden a vigas y 
se asigna las propiedades de esta. 

• [nglYi,ELEMYi,Lyi,senoyi,cosenoyi, VCYi, Yyi,Zyi,NIYi,NJYi,CGYi]=Dat 
os_KLy_3D_lrreg(i,nPortX,ELEM Y,X, Y,Z,NI,NJ) 


Es similar al anterior, pero con los pórticos en sentido Y, por lo que se 
omite su descripción. 

• [Derv]=deriva2(YY,qtiX,qtiY) 


Este programa permite determinar las derivas inelásticas de la estructura 
en sentido XZ y YZ, además genera una gráfica con desplazamientos inelásticos y 
derivas en ambos sentidos donde ei color azul representa ei sentido XZ y rojo YZ 
Los datos son: YY es un vector que contiene las coordenadas en elevación de 
cada nivel; qtiX es un vector con los desplazamientos inelásticos en sentido XZ de 
cada nivel y qtiY es lo propio en sentido YZ. 


4. ANÁLISIS SÍSMICO DE ESTRUCTURA EN 2 D 

En esta sección se presenta ei cálculo de la matriz de rigidez, en sentido 
longitudinal de una estructura irregular de 5 pisos de hormigón armado, en 
muchas estructuras la geometria de los pórticos e incluso las secciones y 
materiales de los elementos de los mismos son iguales para todos los pórticos 
que conforman la estructura espacial y es por ello que resulta interesante solo 
calcular la matriz de rigidez en el plano, es decir de un pórtico tipo de la 
estructura y en base a ello determinar la matriz de rigidez lateral en sentido de 
análisis. La vista en 3D se presenta en la figura 1 a y la vista en planta en la 
figura 1 b donde se indica el pórtico tipo en análisis, tiene 5 pórticos en sentido X 
y 8 pórticos en sentido transversal. 
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Figura 1 a) Vista tridimensional de la edificación; b) Vista en planta de la 
estructura, considerando el pórtico en análisis. 


Se quiere hacer notar que en sentido transversal existen 2 tipos de 
pórticos donde cada uno de ellos son regulares en el plano de análisis; para 
determinar la matriz de rigidez de cada uno de ellos se puede emplear las rutinas 
anteriores de CEINCI-LAB correspondientes a este tipo de planos que son 
detalladas en el artículo “INTERFASE ICEINCI-LAB Y NUEVAS FUNCIONES”. 

Como se menciona anteriormente este ejemplo corresponde a una 
estructura espacial conformada por un pórtico tipo en sentido longitudinal donde 
las secciones y materiales de los elementos dei pórtico indicado en rojo en la 
figura 1 b es constante para la estructura. En la figura 2 se presenta la geometria 
de dicho pórtico y se indica las secciones de cada elemento (columna y viga). 
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Figura 2 Pórtico tipo en sentido longitudinal dei Bloque Estructural; las columnas 
son de 60x60 y las vigas son de 35x60 en los primeros 3 niveles, para niveles 
superiores las columnas son de 40x40 y las vigas de 30x45. 




































68 


Brian Cagua, Julia Pilatasig y Roberto Aguiar 


4.1 Ingreso de datos de geometria de estructura 

La estructura tridimensional con las irregularidades indicadas en la figura 
1 a, puede ser generada desde 3 bloques o partes que componen la estructura 
como se indica en la figura 3 pero dependerá de cada usuário la forma de 
proceder únicamente se debe respetar que los bloques sean generados desde la 
base hacia la cubierta para que sean útiles las rutinas dei análisis de estructuras 
regulares detalladas en artículos publicados anteriormente por los autores. 
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Figura 3 Composición dei Bloque estructural para el análisis sísmico. 

A continuación, se detalla el programa que realiza el análisis sísmico e 
inmediatamente se van comentando la entrada de datos. 
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Figura 4 Parte 1 dei dei Pórtico y código de los datos requeridos. 


Los pórticos de la parte 1 dei Bloque Estructural en sentido longitudinal 
tienen 7 vanos de 6.0 m, 7.0 m y 8.00 m (figura 4); esa información se da en la 
variable “svl”; la altura de cada piso es de 4.0 m para el primer piso y 3.60 para 
el nivel 2, esto se da en la variable “spl”. Los valores de Xol y Yol son ceros 
debido a que este bloque estructural no será desplazado respecto al eje de 
coordenadas X=0, Y=0. 
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Para las partes 2 y 3 que componen al pórtico tipo como se indico en la 
figura 3 se debe ingresar la separación de vanos y alturas de pisos en los 
vectores sp2, sp3, sv2 y sv3 para cada parte, como se detalla a continuación en 
las figuras 5 a y b, esto se lo hace como si fueran pórticos regulares y 
posteriormente para ubicarlos en las posiciones que corresponden a cada pórtico 
se debe ingresar la información de desface en Xo2, Xo3, Yo2, Yo3. Para la parte 
2 se puede notar que Xo2 es cero debido a que el pórtico únicamente debe ser 
desfasado en sentido vertical es decir asignar el valor Yo2 como 7.60 m y para la 
parte 3 que compone el pórtico el desfase debe ser en sentido X e Y por ello se 
tienen los valores de 13 m y 11.20 m correspondientemente. 
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Figura 5 a) Parte 2 dei Pórtico y código de los datos requeridos; b) Parte 3 dei 
Pórtico y código de los datos requeridos. 


Es posible generar el pórtico tipo con más partes de las que se muestra en 
este ejemplo o en algunos casos con estructuras de geometrias más complejas 
será necesario generar más partes que compongan al pórtico tipo y es por ello que 
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se genera el programa Combinar_Vectores con la intención de agrupar la 
información de separación de vanos, alturas de pisos y coordenadas de desfase 
en matrices finales que contengan la información completa de la estructura de 
forma resumida. En la figura 6 se puede observar cómo se emplea el programa 
mencionado donde se combina la información descrita en las matrices svF, spF, 
XoF, YoF. Se puede notar que existen 7 ceros luego de sv3, sp3, Xo3, Yo3, esto 
se debe a que se pueden combinar hasta 10 partes y en este caso únicamente se 
tienen 3, con el objetivo que funcione la rutina de Matlab debe completarse con 
ceros estos “vacíos”. 


•vl, spl , Xol, Yol 1 


Kiun_Cocni»-3, •*-!*■ . • ! il* i 

[ *vF] -Comí»i aat_V*ctor*s , *vl , »v2 . nv3,0 O , O . O , O , O , Q) 

|ipF] s^Zcoioinar Voctorsa |lJun üont>, rpi, splf, ap3 , o, u r 0, ú r J f 3,0) 

(X F|=--mtJÍnjc~V*ct 0, i> 

í YoF] =ç ru r”v*ct or* * ( Nu» ~Coitf> , Yn .,Yn:,yo3,o,o,Q,n,n,o,o> 

Figura 6 Combinar la información de las partes que componen el pórtico tipo. 

En caso de que la estructura sea muy compleja y requiera más de 10 
partes para conformar el pórtico tipo el usuário de CEINCI-LAB puede modificar la 
rutina Combinar_Vectores debido a que el código es abierto y adaptarse a sus 
necesidades particulares de forma sencilla. 

Con los datos ingresados anteriormente y mediante el programa 
geom_2D_lrreg se puede determinar automáticamente la geometria dei pórtico en 
análisis. 


8v2, »p2, Xo2, Yo2 
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%% Geometria de la estructura en 2D 

[nr, nod, np,nELEM,X, Y, NI,NJ]=geom_2D_Irreg(Num_Comb,svF,spF,Xo 
F,YoF); 


4.2 Dibujo de la estructura en 2D. 

Se emplean algunas rutinas típicas para el análisis de grados de libertad y 
otra información relevante para el cálculo de la matriz de rigidez con CEINCI-LAB. 
De igual manera para el dibujo dei pórtico típico se utilizan las rutinas dibujo, 
dibujogdl y dibujoNudElem para obtener los gráficos dei pórtico tipo, incluir los 
grados de libertad dei análisis y enumerar automáticamente nudos y elementos 
con cada rutina (ver figura 7 a, b y c). 
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%% Cálculos de Grados de libertad e información relevante de 
la estructura 

[CG,ngl]=cg_sismo2(nod,nr,Y); 

[L, seno, coseno]=longitud(X,Y,NI,NJ); %long. de @elem, seno y 
coseno 

[VC]=vc(NI,NJ,CG); %Vector de colocación 

%% Dibujo de la Estructura 2D 
dibuj o(X,Y,NI,NJ) 
dibujogdl(X,Y,NI,NJ,CG) 
dibuj oNudElem(X,Y,NI,NJ,CG) 


Cii^Mr*» LnhLf.i.u-il 



< itonaaiH«4 


(b) 
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Figura 7 a) Esquema estructural en 2D; b) Numeración de los grados de libertad, 
con azul se indican las coordenadas principales; c) Numeración de 
Nudos y Elementos, con rojo los nudos y de azul los elementos. 

Al general el pórtico con la recomendación de emplear partes o bloques 
estructurales desde la base hacia la cubierta se puede emplear las rutinas 
anteriores de CEINCI-LAB sin realizar modificaciones adicionales para los grados 
de libertad. En la figura 7c se puede visualizar que los elementos se enumeran 
primero columnas de la parte 1 dei pórtico, luego vigas dei mismo bloque 
estructural, seguido de ello se enumeran columnas y vigas de los bloques 2 y 
luego dei 3, esto es muy importante para asignar las secciones a cada elemento 
dei pórtico, a continuación, se indica como se debe realizar esto. En el programa 

se presenta el valor de la resistência dei hormigón que es de / c ' = 240 / cm 2 y 
la ecuación utilizada para calcular el módulo de elasticidad en T/m 2 . 

%% Propiedades dei material 

fc =240; % La resistência a la compresión dei hormigón f'c 
(kg/cm2) 

E=151000*sqrt(fc); % Calcula el módulo de elasticidad dei 
hotmigón (T/m2) 

%% Secciones 

SECCION=[ 1 0.60 0.60 15 1; 

17 0.35 0.60 27 1; 

45 0.60 0.60 5 1; 

51 0.35 0.60 9 1; 

61 0.40 0.40 7 1; 

69 0.30 0.45 11 1]; 

% SECCION=[i,b,h,nig,ii] 

% i Número dei elemento 
% b Base de la sección dei elemento 
% h Altura de la sección dei elemento 
% nig Número de elementos a generar adicionales 
% ii Incremento en la numeración de los elementos 
[ELEMX]=gelem_portico(SECCIÓN); 
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La matriz SECCION contiene las secciones de columnas y vigas de todos 
los elementos, en la primera fila se escribe 1 0.60 0.60 15 1 puesto que se 
empieza en el elemento 1 con una sección de hormigón armado cuya base es 
0.60 m y altura 0.60 m, a continuación se escribe 15 debido a que los 15 
elementos siguientes también tendrán esa misma sección y el incremento desde 
un elemento anterior a uno siguiente es de 1 en 1, por ello el ultimo valor es 1. 
Para la segunda fila se tiene algo similar con las vigas empezando en el elemento 
17 como se muestra en la figura 7c, luego la base será 0.35m y la altura 0.60m 
generando 27 elementos de la misma sección con intervalos de elementos de 1 en 
1. Para el resto de elemento es algo similar y el usuário puede identificar 
fácilmente el ingreso de secciones en el pórtico de análisis. 

Con la rutina gelem_portico se asigna a cada elemento la base y altura 
correspondiente a la sección de hormigón armado de cada elemento de forma 
automática (en el caso de estructuras de acero se debe asignar área e inércia de 
la sección en lugar de base y altura). 


4.3 Análisis con modelo de un grado de libertad por planta 


Se calcula la matriz de rigidez de cada elemento y se determina la matriz 
de rigidez total (todos los grados de libertad). Para simplificar el análisis se puede 
condensar la rigidez total en los grados libertad principales, es decir los grados de 


en la figura 8. 
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Figura 8 a) Vista tridimensional de la edificación considerando piso rígido; b) 
Esquema simplificado para el análisis sísmico. 


%% Matriz de Rigidez 

[K]=krigidez(ngl,ELEMX,L, seno, coseno,VC,E); %<-MATRIZ "K" 

HORMIGÓN ARMADO (T/m) 

o, 

o. 

KTT=K; %<- MATRIZ "K" DE TODA LA ESTRUCTURA (T/m) 

na=np; 

kaa=KTT(1:na,1:na);kab=KTT(1:na,na+1:ngl); 
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kba=kab'; kbb=KTT(na+1:ngl,na+1:ngl); 
KL=kaa-kab*inv(kbb)*kba; 


Para el análisis en sentido longitudinal se determina la suma de las 
matrices de rigidez de cada pórtico que para el caso es igual a multiplicar por 5 la 
matriz de rigidez dei pórtico tipo. 


KLX = 


npx 

^ KL® 

i=l 


%% Asumiendo que los 5 pórticos son iguales 
KLX=5*KL 

La matriz de masas traslacionales considerando el piso rígido se tiene en 
base ai siguiente código, para este caso la carga reactiva ai sismo se considera 
como 100% de la carga muerta y 25% de la carga viva. 


%% Masas 

% Matriz de masas (T*seg2/m) 

% Se considera W=1*D+0.25*L 
D=0.700; %T/m2 
L=0.250; %T/m2 

ml=(D+0.25*L)* (24*45)/9.8; m2=ml; 
m3=(D+0.25*L)* (24*33)/9.8; 
m4=(D+0.25*L)*(24*20)/9.8; m5=m4; 
M=mdiag(ml,m2,m3,m4,m5); 


Las matrices de rigidez y masas para el análisis sísmico en sentido 
longitudinal tienen una dimensión de 5x5 debido a que la estructura es de 5 pisos 
como se muestra en la figura 8 b, para que el usuário pueda corroborar sus 
cálculos se presenta los valores resultantes para este ejemplo. 


KL X = 


3.44E + 05 
-2.03E + 05 
3.69E + 04 
-1.40E + 03 
1.77E + 02 


-2.03E + 05 
2.68E +05 
-1.21E + 05 
8.70E +03 
-1.04E +03 


3.69E + 04 
-1.21E + 05 
1.14E + 05 
-2.94E + 04 
5.22E + 03 


-1.40E + 03 
8.70E + 03 
-2.94E + 04 
3.92E + 04 
-1.73E + 04 


1.77E + 02 l 
-1.04E + 03 
5.22E +03 
-1.73E + 04 
1.30E + 04 


M = 


84.03 


84.03 


61.62 


37.35 

37.35-1 


4.4 Propiedades dinâmicas 

Con las matrices de rigidez y masas se puede resolver el problema de 
eigen valores y vectores para determinar los periodos de vibración, los modos o 
formas de vibración y frecuencias naturales dei sistema estructural como se indica 
a continuación. 
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%% Propiedades dinâmicas 
[T,fi,OM]=orden_eig(KLX,M); 


Los resultados de periodos de vibración son: 


T = 


rO.767-1 

0.367 

0.184 

0.128 


L0.077J 


4.5 Método Espectral 

Con el programa espectros_nec15MOD se determina las aceleraciones 
espectrales para cada período de vibración, empleando el espectro de la NEC 
2015 obtenido en la ciudad de Quito en suelo tipo D con factores de irregularidad 
en planta y elevación asumidos de 0.9. Con la didáctica pedagógica de que el 
usuário pueda saber que está haciendo en cada etapa dei cálculo se puede 
programar la determinación de los factores de participación modal. 
Posteriormente desplazamientos_modales_CQC encuentra primero los 
desplazamientos en coordenadas de piso, en cada modo de vibración y luego 
utiliza el critério de combinación modal, denominada: Combinación Cuadrática 
Completa, para encontrar las respuestas definitivas en el sentido de análisis. 

Se tiene el cálculo de los desplazamientos laterales y derivas de piso, 
para el sentido longitudinal (para los 5 primeros grados de libertad que se 
muestran en la figura 8 a). Los resultados se indican en la figura 9 a, para la 
componente el espectro de la NEC 2015, que se muestra en la figura 9 b. 

%% Análisis modal espectral 

b= [ 1; 1; 1; 1; 1 ] ; 

R=6; 

fip=0 .9; fie=0 .9; 

is=4; % Suelo D 

iz=5; % Zona sismica 5 = 0.4g 

ip=2; % Region Sierra 

zeda=0 .05; 

[Ad]=espectro_necl5MOD(R,fip,fie,T,is,iz,ip); 

% Ad Aceleraciones Espectrales para cada período de 

vibración en m/seg2 

% Factores de participación modal 

na=np; 

for i=l:na 

FP(:,i)=abs(((fi(:,i))'*M*b)/((fi(:,i)) ' *M*(fi(:,i)))) ; 

end 

[qte]=desplazamientos_modales_CQC(T,fi,Ad,FP,na,OM,zeda) 
for i=l:np 

qti(i)=R* fip* fie*qte(i); 

end 

YY=unique(Y); 

[DervX]=deriva(YY,qti) 
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DESPLAZAMIENTO DERIVA POR PISO 



0 10 20 30 0 5 1 V5 2 25 

Desplazamlento (cm) Deriva de piso (%) 

(a) 


Espectro Sísmico de Aceleraciones NEC 2015 



Período (seg) 

(b) 

Figura 9 a) Desplazamientos laterales y derivas de piso máximas probables 
encontradas, en sentido longitudinal (sentido X); b) Espectro NEC 2015. 


Como se puede observar en la figura 9a las derivas son superiores al 2% 
que se establece como limite máximo permitido en Ecuador por la NEC 2015, 
además debido al cambio brusco de rigidez dei nivel 3 al 4 se puede visualizar un 
cambio en las derivas de piso aumentando radicalmente estas (más de 2 veces en 
las derivas dei piso 4 respecto al 3). 
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5. ANÁLISIS SÍSMICO DE ESTRUCTURA EN 3 D 

En este apartado se presenta el análisis sísmico, en sentido longitudinal 
de una estructura irregular de 5 pisos de hormigón armado, la vista en 3D se 
presenta en la figura 10 a y la vista en planta en la figura 10 b, tiene 6 pórticos en 
sentido X y 3 pórticos en sentido transversal. Se indican los grados de libertad de 
la primera planta, en el centro de masas. 



(b) 

Figura 10 a) Vista tridimensional de la edificación; b) Vista en planta de la 
estructura, considerando los grados de libertad en el nivel 1. 
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Para entender la geometria de la estructura se presenta la geometria dei 
Pórtico A en la figura 11 a, mientras que en la figura 11 b se muestra la geometria 
dei Pórtico B que es igual al C. 

1 2 3 4 S 6 

3.25 I 3.50 3.50 325 3.50 


7*035x20 
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Figura 11a) Pórtico A dei Bloque Estructural; b) Pórtico B y C; las columnas son 

de 30x30 y las vigas son de 35x20. 

















Nuevas funciones dei sistema CEINCI-LAB para el análisis sísmico 


79 


En el programa se presenta el valor de la resistência dei hormigón que es 
de fé = 210 / cm 2 y la ecuación utilizada para calcular el módulo de elasticidad 

en T/m 2 . Nótese que las dimensiones de las columnas se mantienen como 30/30 
cm y de igual manera las dimensiones de las vigas que son de 35/20 cm, todas, 
en los dos sentidos. Las dimensiones de los elementos para el análisis sísmico 
en sentido X se almacenan en la matriz ELEMX y para el sentido Y se cambian 
las dimensiones de base y altura de columnas, registrando estas nuevas 
secciones en la matriz ELEMY. 


Propxedades dal material 

fc «210/ % La tvsiatoncia a U oomprasión ctol hormigón f c (kg/cia2) 

* ml modulo d# • 

Seccionas % Col(30/30); Vigas(35/20) (cm) 

ELEMX-[ 0.3 0*3; 

0.35 0.2]; t Matriz de 2*2 # En la füí 1 va base altura de columnas 
\ en la fila 2 ae eacribe base y altura de vigas 

ELEMY-[ELEMX (1,2) ELEMX( 1 , 1); ELEMX (2,1) ELEMX (2,2)]; * Cambia en sentia Y 

5.1 Ingreso de datos de geometria de estructura 

Para generar la estructura tridimensional con las irregularidades indicadas 
en la figura 10 a, se tiene varias posibilidades y dependerá de cada usuário la 
forma de proceder. En la figura 12 se indica las partes que componen al bloque 
estructural en análisis. 



Figura 12 Composición dei Bloque estructural para el análisis sísmico. 

A continuación, se indica la primera parte dei programa que realiza el 
análisis sísmico e inmediatamente se van comentando la entrada de datos. 


•*t -fl.it.-4l; * * r, j; * i #r un vr.71 
apl -Í2.trl.7f2.T]9 on vmi 

C»..t<jyi -Í0t3.2$f€.75il0.2&il3.&;n) 


* 0 ; Sol - 0 i 
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4.25m 4 00m 



Figura 13 Parte 1 dei Bloque Estructural; a) Vista tridimensional; b) Vista en 
planta; c) Pórtico tipo en sentido X. 

Los pórticos de la parte 1 dei Bloque Estructural en sentido longitudinal 
tienen 2 vanos de 4.25 m y 4.00 m (figura 13 c); esa información se da en la 
variable “svl”; la altura de cada piso es de 2.9 m para ei primer piso y 2.7 para ei 
nivel 2 y 3, esto se da en la variable “spl” para la estructura de 3 pisos. 

El primer dato de la variable CoordYI es 0, que corresponde ai Pórtico 1; 
ei segundo dato es 3.25 m que es la distancia ai Pórtico 2; ei tercero es 6.75 m 
que es la distancia acumulada ai Pórtico 3, para ei Pórtico 4 ei dato es 10.25 m, 
luego se tiene 13.5 m y finalmente se tiene 17.00, distancia ai Pórtico 6 como se 
muestra en la figura 13 b. 

Con los datos indicados se puede generar la estructura regular espacial 
que se muestra en la figura 13 a, tal y como se indica en ei artículo de Aguiar et 
ai. (2019), este bloque estructural puede ser desfasado dei sistema de referencia 
con las coordenadas X=0, Y=0, Z=0 ubicado en la esquina inferior izquierda, las 
distancias para ei desfase de este bloque son Xol =0; Zol =0. 

Continuando la programación se indica un fragmento para la parte 2 dei 
bloque estructural. 

M ClèM* «1 1 « oea itt»d FAATC 2 

•r-' -ií.TJí »írrr«-T nr r * .-/..t .r* >« r .« 

CoordYi -[J.2S.-é.75r lU.ür li.Sr 1?) r Itnçrwaar sti va=»n* cat í mm ercrtanartaa úm laa perfcisa» aii diawasaon X laattrraal 
Nòi -4.J&Í 2>-2 -0.3r \ llllWIMtll ; . . - klIaM áal aiitaaa C mr. |aatlU)> 


Para la segunda parte dei bloque estructural se realiza un bloque regular, 
como se indica en la figura 14 a. En este caso ei bloque estructural tiene 1 vano 
de 4.00 m como se visualiza en la figura 14 b, esa información se da en la 
variable “sv2”; posteriormente se da la altura de piso de 2.70 m que se registra en 
la variable “sp2” y los pórticos en sentido X que componen a esta parte de la 
estructura se indican en ei vector CoordY2; se debe indicar ei desfase dei bloque 
estructural respecto al sistema de referencia con coordenadas X=0, Y=0, Z=0. 
Las distancias de desfase en sentido X se registra como Xo2 y para la altura en la 
cual se debe ubicar a dicho bloque estructural se registra la variable Zo2. 
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(a) 



(b) 


(c) 


Figura 14 Parte 2 dei Bloque Estructural; a) Vista tridimensional; b) Vista en 
planta; c) Pórtico tipo en sentido X. 

De manera similar a los casos anteriores se debe generar la información 
de la parte 3 y 4 dei bloque estructural como se puede apreciar en las figuras 15 y 
16, se debe registrar los desfases correspondientes para generar la estructura 
irregular deseada como se indica en la figura 10 a. 


El código que ejemplifica a estos bloques estructurales se muestra a 
continuación. 


II tláBtá Alt* cot» trrc^uUxida^l 3 

svJ -[C. JÍ 2 i4] | llnaias-i Vôcttür .n li 
=pj \:n»ji«rjc un v tox ccji l- a. 

} *[0a3.25}; i Ir» jceiiar un f* 

Xo) - 0 ; loi - 0 . 3 ;* - i*n - n: 

II Planta a] ta eon t r r»çu 1«ri d»d FAPTX 4 
•V4 *[4) / m v* t, : i f j j 

ap4 •12.1}i tut Oún U « 

. 

X*4 -4.25; ZúC *11; .*.U i-ra 



m 


(b) 


4 2i>m 4.00m 


2.7 m 


rn 


Zo3=8.3/n 


(c) 


Figura 15 Parte 3 dei Bloque Estructural; a) Vista tridimensional; b) Vista en 
planta; c) Pórtico tipo en sentido X. 
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Figura 16 Parte 4 dei Bloque Estructural; a) Vista tridimensional; b) Vista en 
planta; c) Pórtico tipo en sentido X. 

Con la información necesaria de cada parte que compone ai bloque 
estructural se procede a almacenar en matrices los datos referentes a separación 
de vanos “svF”, alturas de pisos “spF”, coordenadas de los pórticos en sentido X 
“CoordYF” y las distancias de desface para X y Z en matrices denominadas “XoF” 
y “ZoF” de cada bloque estructural. Los vectores para combinar pueden ser de 
diferentes dimensiones, pero ai ingresar en la función de Combinar_Vectores 
descrita anteriormente es posible generar matrices a partir de estos vectores. 

4% Combinar las vactoras da cada bloqua qua forman la astructura 30 
Num_Comb-4; I Hace referencia a la cantidad de bloques que forman la estructura 
[svF]-Combinar Vectores<Nu* Comb, svi,sv2,svJ,9v4,0,0,0,0,0,0 ): 

(spFj -Cc*ftL>inai_Vectc>res <Nua_Cainb, spl,sp2, sp3, sp4,0, 0,0,0, 0,0) ? 

[CoordYF]-Combinar_Vcctorea [Num_Carab, C oordYl ,CoordY2, CoordY3,Cc:tjrdY4, 0,0,0, 0,0, 0) ; 
[XoF]=Co«binar_Vectoras <Niun_Comb,Xol,Xo2, Xo3,Xo4,0, 0,0,0, 0,0); 

[ZoF] =Cof*binar_Vectores <Niur._Coinb, Zol, Zo2, Zo3, Zo4,0, 0,0,0,0,0) t 
Caomatna da la astructura an 3D 

[nPortX, nPortY,X, Y, Z,NI,NJ]-geo«n_3D_Irreg(Nuflr._Coinb, svF, spF,CoordYF, XoF, ZoF) s 

Luego de tener la información referente a toda la estructura almacenada 
en matrices se puede generar la geometria de esta de forma espacial, recordando 
que una estructura compuesta por columnas y vigas es descrita por las 
coordenadas de los nudos X, Y, Z, y los elementos que se generan a partir de los 
nudos identificando a cada elemento por el nudo inicial y final de este NI, NJ. 
Además, es recomendable conocer la cantidad de pórticos en cada sentido para 
posteriormente realizar el análisis tridimensional en base a marcos planos. Para 
ejecutar los cálculos de las variables mencionadas anteriormente se emplea la 
función geom_3D_lrreg. 

5.2 Dibujo de la estructura en 3D, de pórticos planos y vista en planta 

Para graficar la estructura en tres dimensiones se utiliza el programa 
dibujo3D el mismo que se muestra en la figura 17. 
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Dibujo do 1* Cstructura JD 
dibu}o3D(X,Y, z, ni # NJ); 

Esquema Estructural 


10 



Figura 17 Dibujo en 3D, de estructura de 5 pisos que se está analizando 
sísmicamente con 3 grados de libertad por planta. 

En algunos casos la visualización de estructuras tridimensionales es 
complicada y es por ello por lo que se desarrolló la función dibujo_2Di con la cual 
se puede graficar pórticos en sentido XZ o YZ dependiendo de la necesidad dei 
usuário, además reporta una vista en planta de la estructura. 

Los datos necesarios para la función dibujo_2Di son las coordenadas de 
los nudos X, Y, Z, los elementos descritos por nudo inicial y final son NI, NJ y el 
número de elementos nelem. La nueva información se refiere a una matriz 
denominada dibujLpian en la primera fila se debe escribir 0 para pórticos en 
sentido X y a continuación de ello la identificación dei pórtico que se desea 
representar, para este caso pórtico 1, 4 y 6; en la segunda fila el primer valor 
correspondiente a 1 es para el sentido Y, seguido dei número 1 y 3 referentes a 
los pórticos A y C, luego para completar la dimensión de la matriz se escribe el 
valor de 0. 


%% Dibujo de pórticos en 2D y viste en plante de le Estructura 
dibuj_plan- [0146; 

1 1 3 0J; 

% dibuj_plan-[sentido(0 es x, 1 es y)/ identificar el pórtico a dibujar] 

nelem=length (NI) ; i Eiement .>s de HA. 
dibujo_J2Di (X, Y, Z,NI,NJ, nelem, dibuj_plan); 
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Figura 18 Pórticos en 2D, estructura de 5 pisos que se está analizando 
sísmicamente con 3 grados de libertad por planta y Vista en Planta. 


5.3 Análisis con modelo de tres grados de libertad por planta 

En la figura 19 a, se presenta ei modelo de tres grados de libertad por 
planta. Nótese que primero se han numerado todos los grados de libertad en 
sentido X, desde ei primer piso hacia arriba; luego desplazamientos en sentido Y, 
finalmente las rotaciones de piso. 































Nuevas funciones dei sistema CEINCI-LABpara el análisis sísmico 


85 


I 


* 






cfc; 

ãi 


CT 

: 1 


Ü 


á) 


rr»5 


nr»4 




ü 


mB 


m2 


ml 


ü 


U 

4 

à*2 

A 

1 


(b) 


Figura 19 a) Modelo de análisis sísmico considerando tres grados de libertad por 
planta; b) Modelación de las masas en piso rígido. 


Se tiene un lazo para el número total de pórticos en sentido X, y lo 
primero que se hace es de la base de datos que se tienen para todos los pórticos, 
se determinan los datos que permiten evaluar la matriz de rigidez lateral de un 
solo pórtico, esto se lo hace con: Datos_KLx_3D_lrreg, con esta función se 
determina, el número de grados de libertad dei pórtico; la base y la altura se halla 
en el archivo ELEMXi, la longitud de cada uno de los elementos dei pórtico se 
encuentra en Lx7; el seno y coseno que permiten encontrar la matriz de rotación 
de coordenadas locales a globales se hallan en: senoxi; cosenoxi; el vector de 
colocación (grados de libertad dei nudo inicial y final) se hallan en VCXi\ y el 
módulo de elasticidad dei material E. 

Para cada pórtico en análisis se reordenan los elementos, las 
coordenadas de los nudos y las propiedades de cada elemento de la estructura 
que se involucra en el análisis. En la figura 20 a se indica la numeración de 
nudos y elementos, con azul los elementos y de rojo los nudos. En estos pórticos 
la numeración de los elementos no necesariamente es primero columnas y 
posteriormente vigas como se acostumbraba a hacer en ejemplos expuesto 
anteriormente. La numeración de nudos y grados de libertad también podría ser en 
un orden diferente al conocido, pero se mantiene que en cada nivel el grado de 
libertad en sentido horizontal de los nudos dei nivel son correspondientes a las 
coordenadas principales (figura 20 b). 
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Esqutma Ettroctural 


li a s a 



Esquema Estructural 



(b) 


Figura 20 Pórtico en sentido X; a) Numeración de Nudos y Elementos; b) 
Numeración de los grados de libertad, con azul se indican las 
coordenadas principales. 

A continuación, en la figura 21 a, se presenta la numeración de nudos y 
elementos, en la figura 21 b se tiene los grados de libertad para un pórtico en 
sentido Y. 


12 

10 

A 

A 

4 

2 

0 



(a) 


(b) 


Figura 21 Pórtico en sentido Y; a) Numeración de Nudos y Elementos; b) 
Numeración de los grados de libertad, con azul se indican las 
coordenadas principales. 
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Con esta información se determina la matriz de rigidez dei pórtico para los 
grados de libertad indicados en la figura 20 b; luego se realiza la condensación 
estática de la matriz de rigidez mediante la solución de ecuaciones lineales. KL. 
Esta matriz se la almacena en una matriz KLXi. 

M AnAliRli 23 Ó9 c»d* poético «n Mfitldo 'XI* 

ÀltuPoi t-unigu* |Z) t % •ndtcuctuiu 

rpT-í «fwjth ifti tnport) - W • -I» pn - pití rtntr nr n 

KLX-reros|npT,np7)r 

for l-ljtiJvjíUC 

.XUtMXl, f _*i # »*naKl f <n7*pnnxt.Vtnri,)t>:l, X*l,«nei,XJlU.C<XM|-D-it-a* XU -ID I rrvgO. nPurtT.RMia,*, V, Z f l»:.IU> » 

M .nrwwo tv 4*1 r^rucc »n imtidH *W 

« cjxt -1 | f^tl . r 7 t j . NI X i r 1 , UX 1 I 

t rti hu] r,Pi 1 f.lnn IX11 < IrJ, VtXi .N tJti f rrTXJ I 


ILHufur t-umtju» |£k;) j « uiv«i i» a i* j . < • >■ X'- 

IXP"l*!n 7 th lAlTTiPort |-lf * Xúnpj, i : - - r- •rt-.w 5 ri : 
% Caütxlbocion ia •looRnt-i ils liúricq f 

x| -kiigiktes (n:|lxi, kLiKXk r L*i, JNinwxi,CEx»«rwaxi, VUJti,£) r 

r.ü-npj «• jdl po? pl-tnti, ooruiiar ar rt*'iip 

k* 4 -X(i:iu, 1 :q 4 I 4 
kAt-r<L j a* f na + lmglXil ; 
t£>a«kab* j 

aür-íC (fu. * 1: nu 1X1, iu« ♦ 1: «4 IX i | / 

r»-Uat\lcb*/ 

KL-tA4*kab*Tf 

%XL-k»* k 4 t* A :nv jfctbí * tb. 

KLl-x«coâ <i»pT f iifj?) 1 KL1 <linp r lsnp|-EL/ 

ft.X» l (í-11 *np7* * 5 1 *tipT r í | -X1.Í r ’-r i -r . rr . t 

KLX-KLX*BliJ 1 j r.u • . - t- rij J«*í ■ t--»r -• 

fc?*d 


KLXi = 


\KL™1 

KL™ 

KL& 

KLW 

KL (S > 

'-KL (6 >1 


Finalmente, en la variable KLX se van sumando las matrices de rigidez 
lateral de los pórticos en sentido X, de tal manera que: 


KLX = 


npx 

Í>“> 

1=1 


Donde la sumatoria se extiende al número de pórticos en sentido X, 
denominada npx en la fórmula 0 nPortX, en el programa. 


Para que se comprueben los resultados, a continuación, se presenta la 
submatriz de la matriz de rigidez en coordenadas de piso, denominada K xx 


K 


XX 


23119.34 -14675.1 

-14675.1 20359.09 

4320.528 -12446.2 

-562.179 2728.48 

28.89 -137.139 


4320. 528 
-12446.2 
15365.83 
-7138.64 
687. 09 


-562.179 

2728.48 

-7138.64 

6513.562 

-1643.57 


28.89 
-137.139 
687.09 
-1643.57 
1069.38SJ 


En el sentido Y, se realiza lo propio, por lo que solo se presenta el listado 
dei programa sin explicación. Solo se indica que se halla 
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KLYi = 


KL™ 

KL (c) - 


npy 

KL W 


i=l 


II Aaáiitii «n 1L' da £*da pótlico ma lanUda “M'* 

Ait«Pcrrt-tinl»jurtZM 1 Mlv»-!»*? ij<» t ---•»• rjctui r 
np7»l«agth(AltuFort|-1# \ Hlh. i qAjuuj >j puo? da v;truotur«i 
KLY-xvroji (np?. apíf j 
to c i-lmPorty 

[oqiíl,ELEHYí. Lyt , awii- r ‘l ,>j_âAaúyi, VCY1« ¥yl.2yl , Hl Vi, MJYi.OLTíl J -Cat..«_KLy_3D_I i iai| U, nK itJt. £L£WY, X, Y, I, H2,BJ)» 



ItltuPrrt-uruqu<? (cyll; * d- -adi pi? • •*« XI 

Up-iWIWtfl fAltuF^I t> - i i I NúlkaXiJ •- 

[ 3 ?]-kiigider inglTi, ELEMYt, Lyl,*©aoyl, cooonoyi, VCYi, E> i 

na-RJ-.; * * j i : r > j nrt r *-- • . f - 1 —,r 

kaa*-K(l jda, Ima) i 
K*fc-X< i 1 ) r 

jtb«“kab*í 

kbfc-K 2 ItvjlYl, J»*-l itejlYlj ê 

KL-Jua«kab*T# 

trL* Md-ltit* invítLt I *kbj 

KLl-m-oa |»f7 ( afiTt t KL* (l:np r l:üpj-KLl 

lanMi-ircpT-itiTrT, ji-Kiix t. . i j,- r.r df í.;- > p 

KLY-KLItia.il liim -1 ut í d.. i-luMwi 1« Kit:. • .. > 

■M 


Para encontrar la matriz de rigidez en coordenadas de piso se utiliza ei 
programa: matriz_es, también se puede utilizar ei programa matriz_es1. Los 
datos de entrada para ei primero de los programas son los siguientes: 

• [KE,rtet,A]=matriz_es(ntot,iejes,np,r,KLT,RT); 

Este programa se explica en algunos artículos escritos por ei primer autor 
dei artículo, se realiza un resumen para que ei lector recuerde los datos y 
resultados de este programa, notot es ei número total de pórticos de la 
estructura; iejes es ei número de pórticos en ei sentido de análisis; np es ei 
número de pisos; r es la distancia dei centro de masas ai pórtico, se da un solo 
valor por pórtico y sirve solo para estructuras regulares en planta y elevación; 
KLT matriz que contiene todas las matrices de rigidez lateral de los pórticos, 
primero en sentido dei análisis símico que para este caso es X, luego en sentido 
transversal que para ei ejemplo será Y; RT es una matriz con las distancias dei 
centro de masas ai pórtico en cada uno de los pisos, esta variable se utiliza en 
estructuras irregulares. 

A diferencia dei artículo “Nuevas Funciones dei Sistema CEINCI-LAB 
para Análisis Sísmico Espacial” en el cual la estructura espacial era regular, en 
esta ocasión los programas presentados son para estructuras irregulares en 
planta y elevación por lo que se debe indicar ai programa r=0. Cuando se ejecuta 
este programa pregunta si la estructura es regular o irregular; en el caso de 
irregulares va a trabajar solo con la variable RT. 
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El programa reporta: KE que es la matriz de rigidez en coordenadas de 
piso; rtet es una matriz que contiene, las matrices con las distancias dei centro de 
masas a cada uno de los pórticos y en cada piso. Finalmente encuentra la matriz 
de compatibilidad A, entre las coordenadas de piso y las coordenadas laterales 
de los pórticos. Aguiar (2012, 2014). 

Datoa para obtener ia matriz de rigidez de la estructura 3D con 3gdl por plao 

Xcml=3,25/2; Ycml = 17/2; Xcm2=Xcml; Ycm2=Ycml;Xcm3=Xcml/ Ycro3=Ycjnl; 

Xcm4-5. 55; Ycm4-7.34; Xcm5-4.25+4.0/2; YanS-3.25/2; 

\ SENTIDO X 

Rl = l Ycml Ycm2 Ycm3 Ycm4 Ycm5]; > F rr, 1 

R2-[Ycml-3.25 Ycm2-3.25 Ycm3-3.25 Ycm4-3.25 Yon5-3.2S]; i:rt;:c 2 
R3=lYcmi-6.75 Ycm2-6.75 Ycm3-t.75 Ycm4-6.75 YcmS-t.75]; • - rt. 

R4-[Ycml-10.25 Ycm2-10.25 Ycm3-10.25 Yon4-i0.25 Ycm5-10.25]; a Pórtico 4 
R$=[Ycml-13.5 Ycm2-13.5 Ycm3-13.5 Ycm4-13.5 Ycm5-13.5]/ - i* 

R6-[Ycroi 17 Ycra2 17 Ycm3-i7 Yon4 17 Yan5-17J? * Mrtico t 
% SENTIDO Y 

RA-fXcni. m3 Xcm4 XcmS]; 

RB«[Xcml-4.25 Xcm2-4.25 Xcm3-4.2S Xc:m4-4.2S XonS-4.25]/ i D 

RO[Xcmi-9.25 Xcm2-3.25 Xan3-8.25 Xcm4-8.25 Xan5-8.25]i % Pórtico C 

% Colocando Ri como diagonal 

Rl-diag(Rl); \ Pórtico 1 

R2-dlag(R2); v pórtico - 

R3< s diag (R3); * Pórtico 7 

R4-dUg<R4); ’ Pr tico 4 

R5«diag(R5); • Pórtico 5 

R6-diag(R6>; % PútUc 5 

RA=diag<RA); % Pórtico A 

RB-diag (RB) ; ' Pórtico B 

RC=diag<RC); * Pórtico C 

RX=[R1/R2;R3;R4/R5;R6]/ RY= [RA> RB; RC] > 

RTX-(RX;RY ]; \ Colocar primei-, la- matuc-. en X y lu-a Y 

ntot=nPortX+nPortY; 

iejes-nPortX; 

r=0; 

%% Analisa s sísmico en sentido X 

KL-[KLXi;KLYl]; 

[KEX,rtet,AJ=matriz_en<ntot,iejen,npT,r,KL,RTX )t 


En el sentido Y, el análisis es muy similar, por lo que no se detalla una 
explicación y se presenta el listado dei programa como sigue: 

%4 Datos para obtener la matriz de rlgldex de la estructura 3D con 3gdl por plao 

iejes=nPortY/ 

rY»0; 

: . : 

_ 

U Analisas sísmico en sentido Y 
KL-[KLYi;KLXi}; 

[KEY, rtet,A]-matriz_e3(ntot,i*3e3,npT, rY,KL,P.TY) i 


5.4 Propiedades dinâmicas 

Como se manifiesta en otros artículos la filosofia de CEINCI-LAB es que 
el usuário también programe algo, para que de esta manera entienda más lo que 
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está haciendo, por esta razón no se utilizo ningún programa para ei cálculo de la 
matriz de masas en coordenadas de piso. 

En este apartado, se indica la carga permanente D considerada en T/m 2 , 
en cada uno de los pisos; para este análisis no se considera la carga transitória L, 
y ei cálculo de la matriz de masas se detalla en ei siguiente fragmento dei 
programa. 


%% Masas 

% Matriz cte nusias (T*3eg2/m| 

% Se considera W=1 * D-í 0*L 

D-0.700; %T/m2 
D2=0.500; VT/m2 

ml-D* (8.25M7)/9.8; m2-ml;m3-inl? 
m4=D2*(8.25*17-4.26*<17-3.25))/9.8; 
m5-D2 é (4 è 3.25)/9.8; 
a=8.25/b=17/ 

Jl-tnl/12* (a A 2*b A 2) ; J2-J1; J3-J1; 
a=4/b=17; 

J4-*m4/12 è (a A 24b A 2); 
a=4;b=3.25; 

J5*nn5/12* (a A 2+b A 2); 

m-mdiag (ml, m2, m3, m4, m5) / J=rodiag (Jl, J2, J3, J4, J5) ; 
M~mdiag(m # m, J); 


El cálculo de los periodos de vibración de la estructura, los modos de 
vibración y las frecuencias naturales dei sistema con análisis en sentido X se lo 
realiza con el programa orden_eig. 


[TX,f lX,aMX]-or<ten_eig<KEX,M> t 

disp<’ Matriz de Piiidez Espacial de la ertructura sentido X ccn *gdl por piso 9 ) 
KEX * Katr J r d® r;gld« * d* 1 Ja '>?tnictura «*n 3P c«.*n 3cjdl por pl.-n 

disp< * * ) 

disp( # Los periodos de vibración en sentido X sor.:') 

TX % aeg 


En la tabla 1 se detalla los resultados de los periodos y modos de 
vibración de la estructura cuando se realiza el análisis símico en sentido X. 

Tabla 1 Propiedades dinâmicas para análisis sísmico en sentido X 



Modo 

Modo 

Modo 

Modo 

Modo 

Modo 

Modo 

Modo 

Modo 

Modo 

Modo 

Modo 

Modo 

Modo 

Modo 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

Período 

(seg) 

0,92 

0,81 

0,71 

0,29 

0,26 

0,23 

0,17 

0,16 

0,14 

0,12 

0,10 

0,09 

0,09 

0,09 

0,07 


- 0,05 

0,01 

0,02 

- 0,12 

0,00 

- 0,09 

- 0,09 

- 0,12 

0,12 

0,03 

- 0,14 

0,13 

- 0,01 

0,00 

0,01 


- 0,14 

0,02 

0,04 

- 0,17 

0,00 

- 0,08 

- 0,03 

- 0,02 

- 0,02 

- 0,02 

0,12 

- 0,16 

0,02 

0,01 

- 0,02 


- 0,20 

0,02 

0,06 

- 0,02 

0,00 

0,05 

0,09 

0,12 

- 0,10 

- 0,01 

- 0,04 

0,13 

- 0,02 

- 0,01 

0,03 


- 0,26 

0,01 

0,01 

0,23 

0,00 

0,13 

- 0,02 

- 0,07 

0,23 

0,06 

- 0,08 

- 0,17 

0,02 

0,03 

- 0,08 


- 0,33 

- 0,05 

- 0,28 

0,58 

- 0,02 

- 0,45 

- 0,40 

- 0,40 

- 0,40 

- 0,16 

0,38 

0,28 

- 0,03 

- 0,10 

0,18 

LJj 2 

0,00 

0,06 

- 0,02 

- 0,01 

- 0,14 

0,02 

- 0,14 

0,09 

- 0,02 

0,16 

0,03 

- 0,01 

- 0,15 

- 0,01 

- 0,01 

Q >0 

0,00 

0,14 

- 0,05 

- 0,01 

- 0,18 

0,02 

- 0,05 

0,02 

0,00 

- 0,09 

- 0,03 

0,01 

0,19 

0,01 

0,02 

1 2 

O co 

§ > 

0,01 

0,20 

- 0,07 

0,00 

- 0,01 

- 0,01 

0,14 

- 0,09 

0,01 

- 0,06 

0,01 

- 0,01 

- 0,15 

- 0,01 

- 0,02 

0,02 

0,25 

- 0,01 

0,02 

0,25 

- 0,03 

- 0,02 

0,01 

- 0,05 

0,27 

0,02 

0,01 

0,18 

0,02 

0,04 


0,03 

0,28 

0,02 

0,03 

0,55 

0,03 

- 0,61 

0,51 

0,17 

- 0,63 

- 0,08 

- 0,03 

- 0,23 

- 0,01 

- 0,05 


0,00 

0,00 

- 0,01 

- 0,01 

0,00 

0,02 

- 0,01 

- 0,02 

- 0,02 

0,00 

- 0,02 

- 0,02 

- 0,01 

0,02 

0,03 


0,00 

- 0,01 

- 0,02 

- 0,01 

0,01 

0,03 

- 0,01 

- 0,01 

- 0,01 

0,00 

0,01 

0,01 

0,00 

- 0,02 

- 0,03 


- 0,01 

- 0,01 

- 0,04 

- 0,01 

0,00 

0,00 

0,01 

0,01 

0,02 

0,00 

0,01 

0,00 

0,00 

0,00 

0,03 


- 0,01 

- 0,01 

- 0,04 

0,00 

0,00 

- 0,05 

0,01 

0,03 

- 0,02 

0,01 

- 0,04 

- 0,02 

0,00 

0,02 

- 0,04 


- 0,01 

- 0,01 

- 0,05 

0,01 

0,00 

- 0,08 

0,02 

0,05 

- 0,07 

0,02 

- 0,26 

- 0,24 

0,03 

- 0,71 

0,19 
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Para el sentido Y se realiza las siguientes operaciones: 

|TY, f tY,OMY|"or<l«n_*lijtKEY,N( i 

dl*p< 'Naltl . Hi~i\ iJi>: r.Mj ... UI !•' 

KEY i M«trlí do ttqldo- d* 1 j «stru 
<U»pC *) 

di»p('t .10 poili l ri 3. vlbt... i:!.n *n 
TY » »*<3 

5.5 Método Espectral 

El usuário debe programar el cálculo de los factores de participación 
modal. Luego de ello llama a dos programas el uno denominado 
espectros_nec15MOD que determina las aceleraciones espectrales para cada 
período de vibración, empleando el espectro de la NEC 2015 obtenido en la 
ciudad de Quito en suelo tipo D con factores de irregularidad en planta y 
elevación asumidos de 0.9. El otro programa que se emplea es el denominado 
desplazamientos_modales_CQC que encuentra primero los desplazamientos 
en coordenadas de piso, en cada modo de vibración y luego utiliza el critério de 
combinación modal, denominada: Combinación Cuadrática Completa, para 
encontrar las respuestas definitivas en el sentido de análisis. 



%% Análisis Modal espoctro de la NEC 2015 
R«4; 

fip—0.9;fie-0.9; 

is-4; % 5ur » D 

12=5; * Zona aiamica 5 = 0.4-j 

ip=2; * Region Sierra 

[Ad]=e3pectro_necl5M0D(R,fip,fie,TX,ia # iz,ip); 

% Ad Aceleraciones Espectrales para cada period£ de vibración en n 
\ Factores de participación modal 

bx«{l;l;1;1;1j;by-[0;0;0;0;0];b 2 =[0;0;0;0;0); 
b=(bx/by/bz); 
na=npT*3; 
for i-l:na 

gamaX(:,l)-abs<< > * *M*b) / < < f IX (: # i)) * *MMf IX (:, i)))) ; 

end 

zeda=0.05; 

( qteX] =desplazamiento3_mcKÍale3_CÔC <TX, f iX, Ad r gamaX, na, 0MX, zeda) ; 
for i-l:npT 

qtiX(i)-R é fip*fie^qteX <i); 

*nd 

YY=AltuPort; 

(DervX)=deriva (YY,«gtiX) 


En la tabla 2 se presenta el factor de participación modal r en los modos 
de vibración y las respectivas aceleraciones espectrales Ad para el sismo con el 
espectro de la NEC15. Se puede observar que el modo 1 es el que tiene 40% de 
participación en el análisis sísmico dejando un 5% para el modo rotacional. 


Tabla 2 Factores de participación modal y aceleraciones espectrales para sismo 

en sentido X 



Modo 

Modo 

Modo 

Modo 

Modo 

Modo 

Modo 

Modo 

Modo 

Modo 

Modo 

Modo 

Modo 

Modo 

Modo 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

r 

5,24 

0,47 

0,99 

1,82 

0,02 

0,89 

0,64 

0,82 

0,66 

0,16 

0,60 

0,52 

0,06 

0,02 

0,05 

r (%) 

40% 

4 % 

8 % 

14 % 

0 % 

7% 

5 % 

6 % 

5 % 

1 % 

5% 

4 % 

0 % 

0 % 

0 % 

Ad 

(m/s2) 

2,72 

3,10 

3,54 

3,60 

3,60 

3,60 

3,60 

3,60 

3,60 

3,42 

3,20 

3,02 

3,00 

2,91 

2,70 


























92 


Brian Cagua, Julia Pilatasig y Roberto Aguiar 


Al final se tiene ei cálculo de los desplazamientos laterales y derivas de 
piso, para el sentido longitudinal (para los 5 primeros grados de libertad que se 
muestran en la figura 19 a). Los resultados se indican en la figura 22 a, para la 
componente el espectro de la NEC 2015, que se muestra en la figura 22 b. 


DESPLAZAMIENTO DERIVA POR PISO 



Desplazamiento (cm) Deriva de piso (%) 

(a) 


Espectro Sísmico de Aceleraciones NEC 2015 



Perodo (seg) 


(b) 

Figura 22 a) Desplazamientos laterales y derivas de piso máximas probables 
encontradas, en sentido longitudinal (sentido X); b) Espectro NEC 2015. 
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Se realiza el análisis modal espectral en sentido transversal, el 
procedimiento es muy similar al realizado en sentido X y por ello solo se mostrará 
únicamente la secuencia de programación sin mayor explicación. 


%% Aiiàlxsxs Modal, espectro de la NEC 2015 
R-4; 

fip-Q.9; fie-0. 9; 
is-1? ~. 1 

iz-5; Zona Llzmi-z* 5 - «.C: 

ip-2; ' Reg: -n fie:ia 

[Adj -esp-ectro_necl5MOD(R, fip, fie,TY r 13 , lz, ip) ; 

% Ad Aceleraeiones Eapectrales para cada período de vibración er 
% Factores de par ticipacion modal 

by-[1;1;1;1;1];bx-[0;0;0;0;0];bz-f 0;0;0;0;0]; 
b-[by;bx;bz]; 
na-npT*3; 
for 1-1: na 

gajna(:,i)-ab3U(fiY(:,i) )' 4 M*b) / < < f 1Y (: , i) ) » 1Y <:, i))) > i 

end 

zeda-0*05; 

(qteY] -desplazamient03_modalea_CQC (TY, fiY, Ad,gama, na # OMY, zeda) / 
for i-i:npT 

qtiY(1)-R*fip*fie*qteY(i); 

end 

YY-AltuPort; 

[DervY]-deriva(YY, qtiY) 


El resultado más relevante tiene que ver con los desplazamientos y 
derivas por esa razón se mostrara en la figura 23 estos resultados dei análisis en 
sentido Y. 


DESPLAZAMIENTO DERIVA POR PISO 



Desplazamiento (cm) Deriva de piso (%) 


Figura 23 Desplazamientos laterales y derivas de piso máximas probables. 
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En la presentación de informes puede ser conveniente realizar una sola 
grafica que contenga los desplazamientos y derivas máximas de piso 
encontradas en ei análisis sísmico en sentido X y Y por ello se presenta la función 
deriva2 cuyo resultado es la figura 24. 



Figura 24 Desplazamientos laterales y derivas de piso máximas probables. 


6. COMENTÁRIOS Y CONCLUSIONES 

Los programas de CEINCI-LAB tienen una capacidad de análisis muy 
poderosa que debe ser aprovechada, en especial porque ei profesional realiza 
conscientemente los cálculos. 

Las estructuras en análisis son muy flexibles, las derivas de piso son 
superiores ai limite máximo de la Norma (2%) y ante un evento sísmico estas 
estructuras son vulnerables estructuralmente y llegarán a sufrir mucho dano. 

La estructura espacial presenta deficiências en ei comportamiento 
estructural que no han sido presentadas pero que evidencian la necesidad de 
realizar un reforzamiento para cumplir como mínimo la normativa y aumentar la 
seguridad de los ocupantes. 

La filosofia dei sistema de computación CEINCI-LAB de ser un programa 
académico se mantiene y con estas nuevas funciones se intenta motivar ai 
usuário para que aprenda a realizar ei análisis sísmico espacial de estructuras 
irregulares paso a paso, se han presentado nuevas funciones con las cuales se 
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dibuja muy fácilmente en 2 D y 3 D, la estructura y la entrada de datos de las 
secciones de sus elementos lo realiza en una forma muy fácil. 

Los programas dei sistema de computación CEINCI-LAB que aqui se han 
presentado se encuentran en Aguiar (2020, 1, 2) y en el siguiente link: 

http://departamentos.uleam.edu.ec/ingenieria-sismica/interfase-en-matlab- 

basica/ 
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